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ABSTRACT: In this paper we continue our program to build a model for high energy soft interactions, 
that is based on the CGC/saturation approach. The main result of this paper is that we have discovered a 
mechanism that leads to large long range rapidity correlations, and results in large values of the correlation 
function R(yi,y 2 ) > 1, which is independent of y± and y 2 - Such behaviour of the correlation function, 
provides strong support for the idea, that at high energies the system of partons that is produced, is not 
only dense, but also has strong attractive forces acting between the partons. 
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1 Introduction 

The large body of experimental data on soft interactions at high energy [1-9], presently, cannot be compre¬ 
hended in terms of theoretical high energy QCD (see [10] for the review). 

In this paper we continue our effort [11-13] to comprehend such interactions, by constructing a model 
that incorporates the advantages of two theoretical approaches to high energy QCD. 

The first one is the CGC/saturation approach [14-20], which provides a clear picture of multi particle 
production at high energy, that proceeds in two stages. The first stage is the production of a mini-jet with 
the typical transverse momentum Q s . Where Q s the saturation scale, is much larger than the soft scale. 
This stage is under full theoretical control. The second stage is when the mini-jet decays into hadrons, which 
we have to treat phenomenologically, using data from hard processes. Such an approach leads to a good 
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description of the experimental data on inclusive production, both for hadron-hadron, hadron-nucleus and 
nucleus-nucleus collisions, and the observation of some regularities in the data, such as geometric scaling 
behaviour [21-25]. The shortcoming of this approach is, the fact that it is disconnected from diffractive 
physics. 

On the other hand, there exists a different approach to high energy QCD: the BFKL Pomeron[26] and 
its interactions[14, 28-35], which is suitable to describe diffractive physics. The BFKL Pomeron calculus 
turns out to be close to the the old Reggeon theory [36], so for calculating the inclusive characteristics of 
multiparticle production, we can apply the Mueller diagram technique [37]. The relation between these two 
approaches has not yet been established, but they are equivalent [34] for the rapidities (In (s/sq)), such that 


Y < 


2 

“BFKL 


In 


1 


A 2 


( 1 . 1 ) 


where A BFKL denotes the intercept of the BFKL Pomeron. As we have discussed [11], the parameters of 
our model are such that for Y < 36, we can trust our approach, based on the BFKL Pomeron calculus. 

This paper is the next step in our program to build a model for high energy soft scattering, based on 
an analytical calculation, without using a Monte Carlo simulation. We discuss the correlation function: 


R(yi,V2;Y) = 


1 


d 2 i7 


1 da 1 da 


&NSD dyidi/2 / \<Jnsd dyi a^sD dY 2 


- 1 


( 1 . 2 ) 


where Y denotes the total rapidity (Y = ln(s/so) and s = W 2 , W is the energy in c.m.f.) and y\ and y 2 
are the rapidities of the produced hadrons. a^sD = &tot ~ &el ~ &sd ~ a dd where a to t(^eh ^sdi^dd) are 
total, elastic, single and double diffraction cross sections. 

The paper is organized as follows. In the following section we discuss the main features of our approach, 
concentrating on the description of diffractive processes. In section 3, we derive the main formulae for the 
correlation functions in our approach, while in section 4 we compare our predictions with the available 
experimental data. 


2 Our model: generalities, elastic amplitude and inclusive production 

In this section we briefly review our model, which successfully describes diffractive[ll, 12] and inclusive 
cross sections[13]. The main ingredient of our model is the BFKL Pomeron Green’s function, which we 
determined using the CGC/saturation approach[ll, 35]. We determined this function from the solution of 
the non-linear Balitsky-Kovchegov equation[17, 19], using the MPSI approximation[38] to sum enhanced 
diagrams shown in Fig. 1-a. It has the following form: 

jessed (T) = fl 2 (1 _ g Xp (_ T ) ) + 2fl ( 1 _ a ) JL- + (! - o) 2 G (T ) 

with G(T) = 1 — ^exp T 0 (2.1) 
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2 

T(s,b) = fa 5(6,m)e a63Aln(s/so) with S(b,m) = ^ e ~ mb 

2n 

In these formulae we take a = 0.65, this value was chosen, so as to obtain the analytical form for the 
solution of the BK equation. Parameters A and fa, can be estimated in the leading order of QCD, but due 
to large next-to-leading order corrections, we treat them as parameters of the fit. m is a non-perturbative 
parameter, which characterize the large impact parameter behavior of the saturation momentum, as well 
as the typical sizes of dipoles that take part in the interactions. The value of m = 5.25 GeV in our model, 
justifies our main assumption, that BFKL Pomeron calculus based on a perturbative QCD approach, is 
able to describe soft physics, since m ftsofti where ft so ft denotes the natural scale for soft processes 
fa so ft ~ A qcd and/or pion mass). 

Unfortunately, since the confinement problem is far from being solved, we assume a phenomenological 
approach for the structure of the colliding hadron. We use a two channel model, which allows us to calculate 
the diffractive production in the region of small masses. In this model, we replace the rich structure of 
the diffractively produced states, by a single state with the wave function ip£>, a la Good-Walker [39]. The 
observed physical hadronic and diffractive states are written in the form 

fa = a \ki + f3 T 2 ; 'tpn = -/? Ti + a T 2 ; where a 2 + fa = 1; (2.3) 

Functions fa and fa form a complete set of orthogonal functions {fa} which diagonalize the interaction 
matrix T 

^ —< ipi 'i/’fc|Tl'i/j/ ip k i >— Ai k Sty $k,k '• (2-4) 

The unitarity constraints take the form 

2ImA itk (s,b) = \A ijk (s,b)\ 2 + G™ k (s,b), (2.5) 

where G"l denotes the contribution of all non diffractive inelastic processes, i.e. it is the summed probability 
for these final states to be produced in the scattering of a state i off a state k. In Eq. (2.5) fas = W denotes 
the energy of the colliding hadrons, and b the impact parameter. A simple solution to Eq. (2.5) at high 
energies, has the eikonal form with an arbitrary opacity 0^, where the real part of the amplitude is much 
smaller than the imaginary part. 

Ai t k(s, b)=i{ l- exp {-fafas, b ))), (2.6) 

Gfkfafa) = 1 - exp (-2 fa fas,b)). (2.7) 

Eq. (2.7) implies that Pf k = exp (—2 Dj ; fc(s, b )), is the probability that the initial projectiles (i, k ) reach the 
final state interaction unchanged, regardless of the initial state re-scatterings. 

Note, that there is no factor 1/2, its absence stems from our definition of the dressed Pomeron. 

In the eikonal approximation we replace Dj fc(s,6) by 

fa k (s,b) = I d 2 b'd 2 b" 9i (?) G“ (T (s,6")) g k (b-b { - ?') (2.8) 
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Figure 1 . Fig. 1-a shows the set of the diagrams in BFKL Pomeron calculus. The wavy double lines denote the 
resulting (dressed) Green’s function of the Pomeron in the framework of high energy QCD, while the single wavy 
lines describe the BFKL Pomerons. In Fig. 1-b we show the net diagrams that include the interaction of BFKL 
Pomerons with colliding hadrons. 


model 

A 

$0 (GeV~ 2 ) 

9i (GeV- 1 ) 

92 (GeV- 1 ) 

m(GeV) 

m\(GeV) 

rri2(GeV) 

P 

ajpjp 

2 channel 

0.38 

0.0019 

110.2 

11.2 

5.25 

0.92 

1.9 

0.58 

0.21 


Table 1 . Fitted parameters of the model. The values are taken from Ref. [12], 


We propose a more general approach, which takes into account new small parameters, that come from the 
fit to the experimental data (see Table 1 and Fig. 1): 

G 31 P / 9i(b = 0) <C 1; m >• mi and m 2 (2.9) 


The second equation in Eq. (2.9) leads to the fact that b" in Eq. (2.8) 
therefore, Eq. (2.8) can be re-written in a simpler form 


is much smaller that b and b' 


^i,k (Sj b) 


^ J d 2 b" G' dressed (T (a,?')) j j d 2 b' gi (?) g k (b - ?) 
(T) j d 2 b' gi (?) gk (b - ?) 


( 2 . 10 ) 


Selecting the diagrams using the first equation in Eq. (2.9), indicates that the main contribution stems 
from the net diagrams shown in Fig. 1-b. The sum of these diagrams[12] leads to the following expression 
for Q iik (s,b) 


n (y ; b) = 


21/ 


d b‘ 


9i[b') 9k[b-b') G 


J l /^dressed 


(T) 


1 + G-sip G diessed (T) 


9i ( b ') + g k (b - b' 


9i (b) = 9i S p (6; rm) 


( 2 . 11 ) 

( 2 . 12 ) 
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where 


S p ( 6 , mi) 


1 „ 

—mi b K i (rriib) 
47T 


^^dressed 


(T) 


J d 2 b G dressed 


(T(s,b)) 


where T (s, 6) is given by Eq. (2.2). 

Note that G dressed (T) does not depend on b and is a function of T = T (s, b = 0) = </>o e°’ 63Ai . 

In all previous formulae, the value of the triple BFKL Pomeron vertex is known: G^jp = 1.29 GeV" 1 
To simplify further discussion, we introduce the notation 


N BK {G^(Y,b)) = a (l-exp(-Gjp(y,&))) + (1 - a) ■, (2.13) 

1 ~r Lrjp [i ,0) 

with a = 0.65 . Eq. (2.13) is an analytical approximation to the numerical solution for the BK equation[35]. 
Gp (Y;b) = gi ( b ) (J dressed (j 1 ). We recall that the BK equation sums the ‘fan’ diagrams shown in Fig. 2. 



Figure 2. A typical example of ‘fan’ diagrams that are summed in Eq. (2.13). 
For the elastic amplitude we have 


a e i{b) — (a 4 A\ } i + 2a 2 f3 2 + /3 4 ^2,2) • 


(2.14) 


To determine the correlation function (given in Eq. (1.2)), we need to know the single inclusive cross 
sections. We have discussed these cross sections in Ref. [13], for the sake of completeness we give the formula 
that describes the Mueller diagram of Fig. 3. 
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9i(b) 

Figure 3. Mueller diagram for the single inclusive cross section. The double wavy lines describe the resulting 
Green’s function of the BFKL Pomerons ( C dressed ). The blobs stand for the vertices which are the same as in Fig. 1. 


+ a 2 /3 2 (/n« Q Y + y'j In® Qy - yj + In™ Q Y + y') In « Q Y - y'j ) 

+ /3 4 In^ Qy + y^j In® Qy - y') \ (2.15) 

where Y denotes the total rapidity of the colliding particles, and y is the rapidity of produced hadron. In W 
is given by 

In® (y) = I d 2 b N bk (g® S (m h b ) G P (y)) (2.16) 

ajpjp is a fitted parameter, that was determined in Ref. [13] (see Table 1). 

3 Two particle correlations 

3.1 Correlations between two parton showers 

The Mueller diagram for the correlations between two parton showers is shown in Fig. 4. Examining this 
diagram, we see that the contribution to the double inclusive cross section, differs from the product of two 
single inclusive cross sections. There are two reasons for this, the first, is that in the expression for the 
double inclusive cross section, we integrate the product of the single inclusive inclusive cross sections, over 
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6, at fixed b. The second, is that the summation over i and k for the product of single inclusive cross 
sections, is for fixed i and k. 

Introducing the following new function, enables us to write the analytical expression: 



a) b) 


Figure 4. The Mueller diagram for the rapidity correlation between two particles produced in two parton showers. 
Fig. 4-a shows the first Mueller diagram, while Fig. 4-b indicates the structure of general diagrams. The double wavy 
lines describe the dressed BFKL Pomerons. The blobs stand for the vertices as shown in the legend. 


Using Eq. (3.1) we can write the double inclusive cross section in the form 

(p(j 2 P ar ton showers 


dyi dy 2 


= / d 2 b < a 4 [ yu b) {y 2 ,b) 


(3.2) 


a 2 /3 2 (yi,b) iW (y 2 ,b) + I^ 1 ) ( yi , b ) I^ 1 ) (y 2 , b)) + /3 4 1 (2 ’ 2) (m, b) I (2 ’ 2) (y 2 , b) 


Comparing Eq. (3.2) with Eq. (2.15) squared, we note the different powers of a and /3 , which reflect 
the different summation over i and k , as well as different integration over b. 


3.2 Correlations in one parton shower: semi-enhanced diagrams 

The main theoretical assumption that we make in calculating the correlation in a one parton shower, is that 
the Mueller diagram technique [37], and the AGK cutting rules[40] are valid. We should note, however, 
that even if the Mueller diagrams provide the correct description of inclusive processes in QCD, the AGK 
cutting rules are not valid for calculations of the correlations in QCD [41, 42]. Nevertheless, we believe 
that, we can neglect the AGK cutting rules violating contributions since, first, they do not lead to long 
range rapidity correlations, which are the main subject of our concern, and second, as we will show below, 
the correlations in one parton shower turn out to be negligibly small. 
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It is instructive to write the expression for the first Mueller diagram in the following form ( see 


,2 one parton shower 2 „y 

asenh — = (ajpjp \n{W/W q)) 2 J2 <x 2 a\ / dY' 


dyi dy 2 


(3.3) 


d 2 b 9i ( b ) G{Y — Y') G 3IP G [y' - l -Y- y^j G (V --Y- y 2 ) 


/ 


2 i/ i\tBK / fc 


do'N' 


9 * (&)G( -Y- yi 


2l/ atBK / „fc 




y*(6)G( -Y-y 2 


The expression for the first Mueller diagram for two parton showers correlation (see Fig. 4-a) has the 
form 

, 2 one parton showers - , 2 2 parton showers 

a U senh _ _ / £,^ a ^ _ _ 

dyi dy 2 J dyi dy 2 d?b 

(fliPiP ln(IT/Wo)) 2 fl a 2 a 2 k f d 2 b ( 9i (b)) 2 G (y'- ±Y 

i=l,k=l d '• 

J d 2 b’N BK (y fe (b) G Qy - yi) j J d 2 b’N BK (^g k (b) G Qy - y 2 ) ) 


(3.4) 


2/1 


G(Y'--Y-y 2 


Comparing Eq. (3.3) with Eq. (3.4) one can see that 


,2 one parton shower 
d ‘°senh 


dyi dy 2 

U(Y,y x ,y 2 ) = 


= T~L (Y. y x , y 2 ) J d 2 b 


rY 


/ dY 

/max(iy-2/i,|Y-j/ 2 ) 


G 3 ip d 2 a ^ P ar t° n showers 
5i ( b ) dyi dy 2 d 2 b 

,G(Y- Y’) G (Y' - \Y - yi) G (Y' — ^Y — y 2 ) 


(3.5) 


G(\Y- yi ) G(±Y-y 2 ) 


H (Y,yi,y 2 ) is proportional to ^Y — yi(yi > y 2 ) in the kinematic region where G(Y) is a constant. At 
small Y it is a constant and is equal to f d 2 bT (Y = 0, b) /A, where A = 0.63 A. Therefore, we can expect 
that the semi-enhanced diagrams can give larger contribution to the double inclusive cross section than the 
production from two parton showers. However, Fig. 5 shows that both the value, and the increase turns 
out to be small in the kinematic region accessible to experiment. Even at ultra high energies, shown in 

Fig. 5-b ,n(Y, yi ,y 2 ) < 0 . 012 . 

On the other hand, the contribution of Eq. (3.5) is small and is proportional to G 3 jp/ 9 i <C 1. Bearing 
in mind that G 3 jp = 1.29 in our approach, one can see that maximum value for max (T~L) ~ 0.012 and the 
values of 

max(%)—~ 1.410 -4 ; max('H)—«1.410 -3 ; (3-6) 

9i (0) 52 (0) 

Therefore, we expect that the contribution of the correlations due semi-enhanced diagrams, is negligibly 
small. 










r H(Y, 0 , 0 ) 




'H( Y, 0 , 0 ) 




Figure 5. Fig. 5-a: The function "H (Y, ?/i, 1/2) versus Y = In (W 2 /Wq ) for J/i = y-i = 0 . The curve 1 shows 
H ( Y, yi, 2 / 2 ) in which all G are replaced by the first order of its expansion with respect to T (Y, b); line 2 is the exact 
calculation. Fig. 5-b is the same as Fig. 5-a but for a large region of Y. In Fig. 5-c H (Y, 2 / 1 , 2 / 1 ) is plotted versus Y 
at different values of y\. H (Y, 0,2/2) is shown in Fig. 5-d as function of I2/2 1 at different energies W. 


The general expression for the double inclusive cross section (see Fig. 6 -b) can be written using two 
new functions S* (2/1,2/2) and S/- (2/1,2/2) defined as 


S* (2/i,2/2) = 


( 3 . 7 ) 


d 2 b’ N bk 5 ( m u b ') G dieaaed Qy + y^j j N BK S (m u b ') G dressed Q 


Y + 2/2 x 


G zip 


Si (221,2/2) = 


( 3 . 8 ) 


J d 2 b' N bk S (mi, b') G dieaaed Qy - ^ N BK ^ 5 (m;, b') G dieaaed Q 


-Y-y. 


It takes the form 


2 1 parton shower 


d~cr. 


semi-enhanced 


dyi dy 2 


= U (Y,2/1,2/2) (a/PJP ln(IY/Wo)) { 2 a 4 S 1 (2/1,2/2) Si (2/1,2/2) 


+a 2 /? 2 (S ' 1 (221,2/2 )-52 (2/1222) + Si (-221,-2/2) S 2 (-221,-2/2)) + 2 / 3 4 S 2 (2/1,2/2) S 2 (2/1,2/2) > ( 3 . 9 ) 
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Figure 6. Mueller diagrams for the rapidity correlation between two particles produced in one parton showers: 
semi-enhanced diagrams. Fig. 6-a shows the first Mueller diagram, while Fig. 6-b gives the structure of general 
diagrams. The double wavy lines describe dressed BFKL Pomerons. The blobs stand for the vertices shown in the 
legend 


3.3 Correlations in one parton shower: enhanced diagrams 


The first Mueller diagram for the correlations from the enhanced diagram is shown in Fig. 7 -a, and has 
the following form 


,2 one parton shower 
d ^enh _ 

dyi dy 2 


f Y 

(ajpjp In (W/W 0 )) 2 ^ a2 i a l / dY ' 


rmin(iy—jn ,\Y-y 2 ) 


(SMfl) 


j d 2 b 9i ( b ) G (Y - Y') G 3IP G (y' - iy - y^j G [y' - l -Y - y^j 


,J d 2 b'g k (b') 


G[-Y + y 2 -Y" 


) g(±Y + Vi -Y”') G 3 jpG(Y") 


(3.11) 


this can be re-written as 


,2 one parton shower 
<P°enh _ 

dyi dy 2 



G 3 ip G 3 jp d 2 a 2 P arton showers 
9 i (b) g k ( b ) dyi dy 2 d 2 b 


(3.12) 


An example of typical diagrams is shown in Fig. 7. The formula, summing all diagrams shown in 
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Fig. 7-b takes the form 


72 1 parton shower 

^ ^"pnhanrpH 


= K (ajpjp In ( W/W 0 )f H (Y, Vl ,y 2 ) U (Y, - yi , -y 2 ) 


dyi dy 2 

|a 4 5 x (yi,y 2 ) S 1 (-yi,-y 2 ) + a 2 fi 2 -S ' 1 (y u y 2 ) S 2 (-y 1} -y 2 ) + S 2 (y u y 2 ) S 1 (-y 1: -y 2 ) 
+ P 4 S 2 (yi,y 2 ) S 2 (-yi, -y 2 ) 


(3.13) 


(3.14) 


where 


K, 



(3.15) 


where Q diessed (Y,b) is determined by Eq. (2.1) and Eq. (2.2). The contributions of enhanced diagrams 
are proportional to the square of the ratios given by Eq. (3.6) and, therefore, they are negligibly small. 



Figure 7. Mueller diagrams for the rapidity correlation between two particles produced in one parton showers: 
enhanced diagrams. Fig. 7-a shows the first Mueller diagram, while Fig. 7-b indicates the structure of general 
diagrams. The double wavy lines describe dressed BFKL Pomerons. The blobs stand for the vertices which are 
shown in the legend 


3.4 Correlation function 


The correlation function R(yi,y 2 ) is defined as 


R(yi,y2) = <?NSD 


d2 2 parton showers j 2 _1 parton shower 1 parton shower 

^ ^ ^ semi-enhanced ^ ^pnhanrpH 


dyi dy 2 


dyi dy 2 


dyi dy 2 


££ - ‘ <-> 
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3.5 Kinematic corrections 


In all our previous equations we assumed that Y = In {W 2 /Wq ) with Wq = lGeV. This assumption 
appears natural for the elastic amplitude, and the cross section of the single inclusive production, but 
it should be re-examined for the correlation function. For this observable, the definition of Y has to be 
modified to account for the fact, that the energy of the parton shower is not equal to W = y/s (s = W 2 ), 
but it is smaller or equal to 5 = W 2 = X 1 X 2 s = X 1 X 2 W 2 (see Fig. 8, where we show the diffractive 
cut of the Mueller diagram of Fig. 4-a). The simplest way to find x± and X 2 is to assume that both 



Figure 8. The Mueller diagram for the rapidity correlation between two particles produced in two parton showers: 
the diffractive cut of the diagram. The double wavy lines describe dressed BFKL Pomerons. The blob stands for 
a p p. 



Pi = P 2 = —Q 2 ^ /tg 0 f t > w here /i go ft is the scale of the soft interactions, A* so ft ~ ^QCD■ In our approach 
the scale of hardness for the BFKL Pomeron Q = m ~^> P so ft- Bearing this in mind, the energy variable 
x\ (X 2 ) for gluon-hadron scattering is equal to 

0 = (x 1 Pi +pi) 2 = -Q 2 + xi2pi ■ Pi; p 2 = -Q 2 - x'i = m P + q 2 - ( 3 - 17 ) 

pi, Pi and X\P\ are the momenta of the gluon, the hadron and the parton (quark or gluon) with which the 
initial gluon interacts. From Eq. (3.17) one has that 

s = xix 2 S = (3.18) 

where M denotes the mass of produced hadron in the diffractive process. For the two channel model, 
it is the mass of the diffractive state. We can use the quark structure function to estimate the typical value 
of x\ = X 2 as it is suggested in Ref. [44]. Using the structure functions at Q 2 ~ 25GeU 2 one finds that 
(|xi|) « 0.3 4 - 0.5. In Fig. 9 the values of R (Y, 0, 0) are plotted for Y = In (x\W 2 /Wq) as function of x\. 

3.6 Correlation in one parton shower: emission from one BFKL Pomeron 

In addition to the sources of correlation that have been discussed above, we need to take into account the 
correlation between two gluons emitted from one BFKL Pomeron (see Fig. 10). At large yi 2 = \yi — 2 / 2 1 
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R(Y, 0 , 0 ) 



- W = 13 TeV 

- W = 7 TeV 

- W = 1.8 TeV 
- CMS data 


Figure 9. R (Y , 0, 0) versus x\ for Y = In {x\W 2 /W$ ) for different energies W. The red line indicates the moment 
C 2 = (|?7 2 1)/ ((|n 2 |)) 2 = 2 for the window in rapidity —0.5 < y < +0.5, measured by the CMS experiment [ 8 ]. 

the diagram of Fig. 10 induces long range correlations in rapidity, however, at small y 12 this emission is 
suppressed, and we do not expect a large contribution from this source. 

The contribution of this diagram can be written in the form 

i? BFKL (yi,y 2 ) = (3.19) 

_ Ei,fc=i a i a k?i - Vi) G dieBBed (y 12 ) T k (\Y + y 2 ) _ 

NSD Ei,fc=1 on a k Ti {\Y - Vl ) T k {\Y + yi ) Ej fc=1 on a k T t (±T - y 2 ) {\Y + y 2 ) 



Figure 10. Mueller diagram for emission of particles from one BFKL Pomeron. 


where 

r i(y) = J d 2 bN BK ( gi (b)G“{y )) 


(3.20) 
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For small values of the argument Tj (y) V - << - > f d 2 bgi ( b ) (5 dressed (y) Eq. (3.19), has no dependence on 
yi2, leading to long range rapidity correlations. However, it turns out that the exact computation, leads 
to very small values of R BFKL (yi,y2): approximately 0.2 -f- 0.4 % of the contributions from the sources 
discussed above. 

3.7 Short range rapidity correlation 

Besides long range rapidity correlations, the emission from one BFKL Pomeron, as well as the hadronization 
in one gluon jet, can lead to short range correlations in rapidity. Unfortunately, at present, this contribution 
cannot be treated on pure theoretical grounds, as it is involves, confinement effects. To estimate this 
contribution, we introduce the Mueller diagram, shown in Fig. 11, where we describe this correlation by 
the phenomenological constant ajpjR, and introduce the correlation length A ~ 2. In the diagram of 
Fig. 11 for the zigzag line we have flp K exp (-^) (yi2 = |yi = y2 1 )• Our estimate for A stems from 
Reggeon phenomenology, in which the zigzag line describes the contribution of the secondary Reggeon, 
with a propagator exp (—(1 — Off?,(0))yi2) and o/r(0) ~ 0.5. 



9j(b) 


Figure 11. Mueller diagram for the short range rapidity correlations. Wavy double lines denote BFKL Pomerons, 
zigzag lines describe the exchange of the secondary Reggeon trajectory. 

The contribution of this diagram takes the form 

^.short range (y i ,y 2 ) = ( 3 . 21 ) 

_ Ei,fc=i Qi r, (ly - yi ) e-^ r fc (iy + V2 ) _ 

GNSD £?,*=! o. k Yi (iy - yi) T k (iy + yr) EL=i ^ ^ {±Y - y 2 ) (ly + y 2 ) 

4 Predictions and comparison with the experiment 

Fig. 9, shows that the correlation function increases with energy, and becomes rather large (of the order of 
1) at W = 7 TeV. This qualitative feature is in agreement with the experimental data from the LHC. The 


/ QjPiR 
\ajpjp 
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first set of data, is the multiplicity distribution measured by the CMS collaboration^]. In particular, the 
value of C *2 = (|n 2 |)/ ((l ra2 |))^ turns out to be very close to 2, for the window in rapidities —0.5 < 17 < 0.5. 
Since C 2 = R(0,0) + l + l/dN/dr]\ ri= o where dN/ dry ]= Q denotes the multiplicity at rj = 0, and at W = 7TeV 
it is equal to 5.8, while R (0,0) = 0.82*. 

The second set of the data is the measurement of the double parton interaction (DPI) [45]. In the LHC 
experiments, the double inclusive cross sections of two pairs of back-to-back jets with momenta px.\ and 
pr, 2 , were measured with rapidities of two pairs (y\ and j/ 2 ) close to each other (y\ ~ 1 / 2 ) ■ These pairs can 
only be produced from two different parton showers. The data were parameterized in the form 

da m da da ., _, 

---5- =- (4-1) 

dyid 2 p T) idy2d 2 p T ,2 2 cr eff dyid 2 p T ,i dy 2 d 2 p T ,2 

where m = 2 for pairs of different jets, and m = 1 for identical pairs. One can calculate the rapidity 
correlation function using Eq. (4.1) 


R(yi,y2,PT,i,PT,2) 


1 da 

&in dyid 2 p T ,i dy2 d' 2 p T ,2 

1 _ da _1_ da 

Pin dyid 2 p T , 1 cr in dy 2 d' 2 p T ,2 


1 = - 1 „ 2 

a eff 


(4.2) 


For the above the estimates we use cr e ff — 12 - 15 mb (see Refs. [45]) and a m = — a e i — a s d — a<id ~ 50 

mb for the energy W = 7TeV (see Ref. [12] and references therein). These data confirm that at high 
energies we are dealing with a system of partons that have a large mutual attraction . The fact that we 
predict a smaller correlation than in this experiment, does not discourage us, since the correlation function 
in Eq. (4.2) differs from the one that we calculate (see Eq. (1.2)). 

The forward-backward correlation has been measured by the ATLAS collaboration in Ref. [9]. The 
observable that was used in Ref. [9] , differs from the correlation function R ( 171 , 112 ), and it can be re-written 
as 


Pfb 


R(m,m) 


(4.3) 


y/R(m,vi) R(v2,V2) 

The value of ~ 0.666 [9] indicates large correlations, but it is difficult to compare with our estimates, 
since ATLAS introduced a specific selection: the px of all produced particles should be larger than 500 MeV, 
while R ( 17 1 , 172 ) is defined as integrated over all momenta. 


Using R( 0,0), we estimate the values for C 2 for W = 13TeU, and using the formula for the negative 
binomial distribution: 


a 


( r V T (n + r) 
\n + r J n! T (r) 


n 


n + r 


n 


(4.4) 


we obtain the multiplicity distribution in the rapidity window —0.5 < 17 < 0.5 shown in Fig. 12. In Eq. (4.4) 
n = dN/dp 1^=0 which was calculated in our paper [13]. From Fig. 12 we expect a violation of the KNO 
scaling behavior[46]. Accordingly to KNO scaling a n /a = F (n/n) with n = J_°q 5 drjdN/dp = dN/dp 1^=0- 


*In this section we use pseudo rapidity y instead of rapidity y, since this variable is used to display data from LHC 
experiments. We recalculate y = h ( y ) where function h is taken from our paper [13]. 
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Figure 12. 


tr n /cr versus multiplicity n/n at different energies W in the rapidity window — 0.5 < y < 0.5. 


It turns out that R ( 2 / 1 , 2 / 2 ) at fixed energy, depends neither on y\ nor on y 2 , giving a perfect example 
of long range correlations in rapidty. To understand why we have these features, it is instructive to start 
from Fig. 4-a at small values of Y. In this kinematic region we can replace G' dressed —>■ G bare = T ( Y ) and 
N BK ( Y,b ) —» gi ( 6 ) G bare . After simple algebra, the correlation function is equal to 

2 


1 


R(yi,U2) = —t vnsd 




d 2 b' gi 


9i 


b-b' 


or 


( K Jd 2 bd 2 b'g 1 (b) gi (b-b’))' 


- 1 


(4.5) 


In Eq. (4.5) we use the fact that in our model g± 3 > <72 • Eq. (4.5) leads to a correlation function that does 
not depend on 2/1 and 2 / 2 - 

On the other hand, at very large Y, N BK ( Y,b ) —>• @(R(Y) — b ), where 0(6) is a step function. 
Plugging in this simple expression, we obtain 


R(yi,V2) 


f d 2 b 0 (R (ly - y 2 ) -b)G(R (jY - yi ) - b) R 2 (iy - yi ) R 2 {\Y - y 2 ) 

( 1 y _ yi ) R 4 (1 Y - y 2 ) 

_ PNSD _ _ ^ 

irR 2 (|y - yi) 


- 1 
(4.6) 


R (y) in Eq. (4.6) denotes a typical impact parameter at large Y, which is proportional to Y t. Recall, 
that at high energies, all components of the wave functions in the two channel model, give the same 
contribution. This is the reason why we do not have an extra factor which depend on a and (3. 


Eq. (4.6) shows the logarithmic dependence on y±. Using Eq. (4.6) we can estimate the 271 ( 2 / 2 ) depen¬ 
dence of R ( 2 / 1 , 2 / 2 )) calculating 


( \b 2 (Y)\) = 


f b 2 d?bN BK [gi (6) G diessed 


f d 2 b N BK [ gi (6) G dressed 


(4.7) 


t\Ve trust that our use of the same notation for the correlation function and typical b, will not confuse the reader 
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Figure 13. (| b 2 (Y) |), defined in Eq. (4.7), versus Y. 


From Fig. 13, in which we plot the results of our calculation, one can see, that only at large Y > 11 does 
(| b 2 (Y) |) start showing visible Y dependence. Two vertical dotted lines mark the widow in rapidity, which 
is essential in the calculation of the correlation function at W = 7TeV for —2 < 2/1 >2/2 < 2. We can 
expect a change of R by 2%. The actual calculation gives even less. Using N BK (Y, b ) —>• 0 (R (Y) — b) 
gives R? (Y) = 2(| b 2 (Y) |). Using this relation, we estimate 7? ( 2 / 1 , 222 ) as 


R{yi,y2) 


J-_ &NSD _ 

a 4 2 tt(|6 2 (\Y - yi ) |) 


1 


(4.8) 


At W = 7 TeV from Eq. (4.8) we find that -R ( 2 / 1 , 2 / 2 ) = 1-647, while the exact calculation give 1.437 
(see Fig. 9). At W = 13 TeV this simple formula leads to ii ( 271 , 2 / 2 ) = 1-72, versus 1.64 from the exact 
calculations (see Fig. 9). 

The correlations, measured by the ATLAS collaboration [9], at first sight contradict both our esti¬ 
mate and the CMS data, regarding the multiplicity distribution. We first check Eq. (4.3). The measured 
observable has the form [9] 

(fa/ “ ( n /» (n b ~ K»> 

Pfb ~ 7 l 4 - y ) 

V (( n / _ <"/)) ) i( n b - (nb)) ) 

The numerator of Eq. (4.9) can be written as R ( 2 / 1 , 2 / 2 ) A 2/1 A 272 , where A yi is the interval of rapidities 
where the hadrons are measured. However, at the same value of rapidity R ( 2 / 1 ,2/i) (4\?/i) 2 corresponds to 
(n(n — 1)) — ( n ) 2 . Therefore, the expression for 


<("/ - (n f )) 2 ) 


( 2 / 1 , 2 / 1 ) + 




(A2/1) 2 
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which leads to the following formula for p 


Pfb 


R(yi,V2) 


R(yi,yi) + 


dy\ 


R(y2,y2 ) + dJV 




(4.10) 


Taking Ayr = 0.5, we see that the first element of the Table 2 is equal to 0.7, which is in good agreement 
with the experimental value 0.666 ± 0.011. 


Forward r] interval 

Backward 77 interval 

0.0 - 0.5 

0.5 - 1.0 

1.0 - 1.5 

1.5 - 2.0 

2.0 - 2.5 

0.0 - 0.5 

0.666 (0.70) 
± 0.011 

0.624(0.643) 

± 0.011 

0.592 (0.599) 
± 0.011 

0.566 (0.565) 
± 0.012 

0.540 (0.539) 
±0.013 

0.5 - 1.0 

0.624 (0.667) 
± 0.011 

0.596(0.618) 

± 0.011 

0.574 (0.580) 
±0.012 

0.553 (0.550) 
±0.013 

0.530 (0.527) 
±0.014 

1.0 - 1.5 

0.594 (0.640) 
±0.011 

0.576(0.596) 

±0.012 

0.560 (0.563) 
±0.013 

0.540 (0.537) 
±0.014 

0.518 (0.516) 
±0.014 

1.5 - 2.0 

0.571(0.615) 

±0.012 

0.557(0.577) 

±0.013 

0.544(0.548) 

±0.014 

0.526 (0.525) 
±0.014 

0.503 (0.508) 
±0.016 

2.0 - 2.5 

0.551 (0.593) 
±0.013 

0.540 (0.560) 
±0.014 

0.527 (0.535) 
±0.014 

0.507(0.515) 

±0.016 

0.487(0.499) 

±0.018 


Table 2 . Multiplicity correlations for events at = 7 TeV for events with a minimum of two charged particles in the 
kinematic interval px > 100 MeV and \r]\ < 2.5 for different combinations of forward and backward pseudorapidity 
interval. The data is take from Ref. [9]. The numbers in parentheses are our estimates. 

To describe our results given in Table 2, we need to also take short range rapidity correlations into 
account. In this table, in parenthesis, we have our estimates, which we obtain on describing the correlation 
function in the form 


R{yi,V2) = 


R 


short range 


(yi 


R UmR ra "ge {y u y 2 ) + ^ shortranse (y 1 ,y 2 ) - i? shortrange (0,0)J; 

, _ _ r, (ly-mp-Tr, {jY + m ) _ 

m “ O (Jr - S1 ) r, (ir + »,) r, (|r - S2 ) r, (iy + V2 ) 


(4.11) 


with a = 0.7 and A = 2. In Eq. (4.11) yi 2 = |yi — 2 / 2 1 an d we restrict ourself by the contribution of the 
state "1" in Eq. (3.21), since g\ y 2 . We assumed that i? short range (0,0)=0 since at yi 2 = 0 Eq. (3.21) 
leads to the long range correlations which we have calculated in section 3.4. One can see that the agreement 
is not perfect, but it demonstrates that the ATLAS data can be reproduced, by including the short range 
correlations. 
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5 Conclusion 


The main result of this paper is that in our model, which is based on the CGC/saturation approach, we 
have discovered a mechanism that produces large, long range rapidity correlations at high energies. The 
large values of the correlation function R(yi,y 2 ) > 1 at high energies, lends strong support to the idea, 
that at high energies the system of partons that is produced, is not only dense, but also has strong attractive 
forces acting between the partons. The resulting long range rapidity correlations are independent of y 
. This prediction is in direct contradiction to the estimates from the soft Pomeron based model that we 
made (see Ref. [43]). In that model the correlations from one parton shower are larger than from the two 
parton showers, and they led to the y\ .2 dependence. Scrutinizing our formulae we found that the main 
reason for the smallness of the correlation in one parton shower that we observed in our approach, stems 
from the most theoretically reliable part of our model: from the expression for the dressed BFKL Pomeron 
Green function. 

We demonstrated that our model is able to describe the LHC data, eminating from the CMS and 
ATLAS collaborations. These data are certainly insufficient for a thorough analysis of the details of our 
approach, but they confirm that the long range rapidity correlations are large at high energies. The 
prediction for W = 13 TeV is shown in Fig. 9. The correlations should increase with the energy and the 
results measurements at W = 13 TeV will clarify the situation. 

In general we believe that this paper is the next natural step, in building a model capable of describing 
soft high energy interactions based on the CGC/saturation approach. 
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